Entanglement in quantum computers described by the XXZ model with defects 
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We investigate how to generate maximally entangled states in systems characterized by the Hamil- 
tonian of the XXZ model with defects. Some proposed quantum computers are described by such 
model. We show how the defects can be used to obtain EPR states and W states when one or two 
excitations are considered. 
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I. INTRODUCTION 

Since qubits are two level systems, they are natu- 
rally modeled by spin-1/2 particles. Understanding spin 
chains is therefore very useful in the study of quan- 
tum computers (QCs). Interaction between qubits corre- 
sponds then to interaction between spins. One of the ma- 
jor problems of condensed matter based quantum com- 
puters is that the interaction between qubits cannot be 
turned on and off when desired, and the quantum com- 
puter eigenstates soon become a linear superposition of 
a large number of noninteracting multi-qubit states 0. 
However, when performing computations, we would like 
to operate with well defined states, in other words, we 
would like to entangle just some specific states. In or- 
der to do so, we recur to an important characteristic of 
most proposed QCs, that the energy difference between 
the qubits states is large compared to the qubit-qubit in- 
teraction and can be individually controlled 0, U • ^ n the 
QC based on electrons on helium, for example, the level 
spacing of each qubit is controlled by electrodes placed 
beneath the helium surface Q- The possibility of indi- 
vidually control qubit energies allows us to entangle just 
some selected qubits by tuning them in resonance. 

A fundamental requirement for the realization of quan- 
tum computation, quantum teleportation and some pro- 
tocols of quantum cryptography is the generation of 
highly entangled quantum states. The maximally bi- 
partite pure-state entanglement is identified with the 
Bell or also called Einstein-Podolsky-Rosen (EPR) state 
(1/V2)(|10> + |01) ) 0. Diir et al showed that there are 
two different kinds of genuine tripartite pure state entan- 
glement: the maximally entangled Greenberger-Hornc- 
Zeilinger (GHZ) state Q and the so called W state 0. 
The W state is the state of three qubits that retains 
a maximal amount of bipartite entanglement when any 
one of the three qubits is traced out. It is written as 
(1/V3)(|100) + |010> + |001)). 

There is a vast list of references dedicated to the prob- 
lem of entanglement. There have been attempts to char- 
acterize qualitatively and quantitatively the entangle- 
ment properties of multiparticle systems 0- Several 



proposals of how to prepare entangled states in differ- 
ent kinds of systems have also been presented . Some 
papers investigated the entanglement between spins in a 
one dimensional Heisenberg chain |9lll0j. which is similar 
to what we intend to do here. But contrary to these last 
cited papers, we consider a chain with defects and study 
their role in entangling states. In a system where all 
qubits are in resonance but one, a defect corresponds to 
the qubit whose level spacing is different from the others. 

In this paper we investigate how to entangle selected 
qubits in a system described by a strongly anisotropic 
one-dimensional XXZ model with defects. This is the 
model used to describe the quantum computer based on 
electrons on helium p|. Since the coupling falls down 
quickly with the interqubit distance, we consider only 
nearest neighbor interaction. The ground state of the 
system corresponds to all spins pointing down and exci- 
tations correspond to spins pointing up. The interaction 
can only move excitations one site to the left or to the 
right, so the number of excitations is constant. We ana- 
lyze systems with one or two spins up. In the case of one 
single excitation we show how the defects of the chain 
can be used to maximally entangle two qubits. In the 
case of two excitations we show how EPR states and W 
states can be obtained. 



II. GENERATION OF MAXIMALLY 
ENTANGLED STATES OF SELECTED QUBITS 

The Hamiltonian of the XXZ model with defects is 
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where h — 1, <r z ' + '~ are Pauli matrices, and e n gives 
the energy difference between the two states of qubit n. 
There are N qubits. Here we consider a spin chain with 
free boundaries, which explains why the second sum runs 
over n — 1, N — 1. In strongly anisotropic systems 
(such as the QC based on electrons on helium), the pa- 
rameter A is much larger than B. The last term in the 
Hamiltonian, H^p, is responsible for the propagation of 
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the excitation. 
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The ground state corresponds to all spins pointing 
down and its energy is Eq = — Yl n =i £ n/2 + (N— l)£?A/4, 
which we will set equal to zero. Moreover, to simplify our 
analysis we consider only positive values for the parame- 
ters of the Hamitonian. 

To address the different states of the system we use 
a notation that is common in the study of spin chains 
with the Bethe ansatz ^J, |^. The state correspond- 
ing to one single excitation on site n, that is | |iJ,2 
• •■ Iri-itnln+i •■■ In), is simply written as 4>(n). The 
state of two excitations, one on site n and the other one 
on site m, is <f>(n,m), which is a simplified notation for 
I li ••• Tn ••• Tm ••• In)- 



A. One excitation 

Let us first examine the case of just one excitation. As- 
sume that there are only two defects, whose level spacings 
are s + g, while the level spacing of all other qubits is e . 
By choosing g much larger than the interaction strength 
B, we generate maximally entangled states correspond- 
ing to linear combinations of the two defects. 

The energy of states that have the excitation on any 
site but on the defects lies within the band £i±B, where 
E\ = Eq — BA 0, If g is much larger than B, 

an excitation on one of the defects will have energy out 
of the band. Therefore the two resonant defects can be 
treated separately from all other states by perturbation 
theory. It is as if we were working with only two sites. 
An excitation initially created on one defect will only 
oscillate between the two defects. All intermediate states 
for the excitation to go from one defect to the other are 
virtual states. The frequency of these oscillations depend 
on the distance between the two different qubits. Their 
separation determines in which order of the perturbation 
theory they are connected. 

The two states corresponding to superpositions of an 
excitation on a defect on site no and an excitation on a 
defect on site toq are the following EPR states 



i>± = [<t>(no) ± <P{™a)], 



(2) 



where </>(n ) = I Tn 4m ) and 4>(m ) = | |„ tm ) 

If the two defects are next to each other (mo = n + 1), 
the energies of these two entangled states in first order 
are 



E± = £ 1 +g±B/2 



(3) 



We just need to diagonalize a two dimensional sub ma- 
trix, whose diagonal elements arc £ i + g and off diagonal 
elements are B/2. After shifting the energy levels by a 
second order correction B 2 /4(g + B/2), they agree very 
well with a complete numerical diagonalization of a long 
chain with g much larger than B. 



If the second defect is located on site no + 2 we have 
to go straight to second order and diagonalize a matrix 
whose diagonal elements are E\ + g + B 2 /(2g) and off 
diagonal elements are B 2 /(Ag). The two energies become 
much closer, 



E+=£i 
E- = £i 



g + 3B 2 /(4g) 
9 + B 2 /(4g). 



(4) 



The more distant the two defects are the closer the ener- 
gies of the two entangled states will be, since we have to 
go to higher orders to find them. 

Suppose an initial state is prepared which has an exci- 
tation on the defect uq. Let us now see how long we have 
to wait for it to become a maximally entangled state such 
as given by Eq.©. This excitation oscillates between the 
defect riQ and the defect mo- The probability to find it 
later in time on site no is 



1 + cos[(£+ - EJ)t] 



while the probability to find it on toq is 



^(m )(*) — 



1 - cos[(£; + - E-)t] 



(5) 



(6) 



It is seen from (JjJJ and JJjJ that the period of oscil- 
lation of the excitation between the defects is inversely 
proportional to the energy difference of the states "0+ and 
ip—. Such period depends on the number /i of sites be- 
tween the defects as T M = T (2g/By, where T = 2-ir/B. 
The maximally entangled states @ are obtained when 
P<t>(n ) = P<f>(m ) = 1/2- The further the defects are from 
each other the longer we will have to wait for an EPR 
type of state to be created. It is clear, however, that 
by tuning two qubits in resonance with energies very dif- 
ferent from all the others, we can entangle even remote 
qubits. 

At the moment where a maximally entangled state is 
created, in order for it to be kept this way, the two defects 
have to be detuned. The difference in energy between 
these two excited qubits should become much larger than 
the interaction strength between them (of course the de- 
fects and the other qubits are still completely out of res- 
onance) . How fast this detuning should be done depends 
on how much close we want to keep our state from a 
perfect EPR state. 

Similarly, a W state can be built with three resonant 
defects, but we postpone the description of how to create 
such states to the next section where we have the more 
interesting case of two excitations. 



B. Two excitations 

In order to create EPR states and W states with two 
excitations we make use of the anisotropy of the system. 
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Because A is much larger than B, two excitations next 
to each other have energy much larger than any state 
where they are separated. The energy of any two free 
excitations lies inside the band 2£i±2B, while the bound 
pair states have energy within a much narrower band 
2£i + BA + B/2A ± B/2A 0, D3- As a consequence, 
the bound pair states can be treated separately from all 
other two-excitation states |l2j. They are connected in 
second order of perturbation theory, which explains the 
narrow bandwidth B 2 /(BA) = B/A. Any intermediate 
and dissociated state is a virtual state. 

Suppose that there is only one defect on site no with 
level spacing e + g. If g is much larger than B/2A, the 
bound pairs with one excitation on the defect have energy 
out of the narrow band. They form the EPR states 



= -^=[0(n o - l,n ) ± 0(no,no + 1)] (7) 

where 0(n o - l,n ) = | T„ -iT„ |„ +i) and </>(n ,n + 
1) = I Ino-iTnoTno+i)- And tneir energies are 



E± = 2£x+g+BA- 



B 2 
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4A 4{BA + g) A(BA + g) 



(8) 



By preparing an initial state with one excitation on 
site no — 1 and the other on the defect site no, following 
(JSJ) and © , we will obtain an EPR state at every instant 
of time 



t k = 2{BA + g)[n/2 + kir]/B 2 



(9) 



where k is an integer number. 

Using the anisotropy, several other types of EPR states 
can be created. With three defects on sites no, no + 1 
and no + 2, if they all have the same level spacing Eq + g 
and g S> 5/(2 A), we would have linear combinations of 
the bound pairs (f>{n 0l n + 1) = | Tnotno+i-Ua+2) and 
0(n o + l,n + 2) = I lnotno+itno+2)- As mentioned 
above, these states are connected in second order of per- 
turbation theory. 

Another EPR state that can be created with those 
three defects involves the states </f>(no,n + 1) = 
) and (/>(n ,n + 2) = | Tno-Lo+i Tno+2/ • 
The level spacing of the defect on site no + 2 is now 
e o + .9 + BA. The difference BA from the other two 
defects allows the entanglement between the bound pair 
</>(no,no + 1) and the state 0(no,no + 2). These states 
are connected in first order in B. The advantage of this 
entanglement is that the period of oscillations between 
an initial state </>(no,no + 1) and the state 0(no,no + 2) 
is much shorter than oscillations between bound pairs. 
On the other hand, to guarantee that </>(no, n a + 1) and 
</>(no, no + 2) are the only two states of the entanglement, 
g has to be larger than B, instead of just larger than 
5/(2 A) as in the bound pair case. 



This anisotropic chain with defects can also be used 
to create a W state. In order to do so we choose four 
equal defects on sites no , no + 1 , no + 2 and no + 3 with 
level spacings £0 + 5 and g much larger than B/(2A). 
The bound pairs on the defects are much higher in en- 
ergy than any other state and they can be treated sep- 
arately. It becomes a good approximation to say that 
three of the eigenfunctions of the total Hamiltonian 
will correspond to linear combinations of the bound pairs 

0(n o ,n o + 1) = I Tn Tno+l4no+24n +3), 0("O + l,n + 
2) = I 4n Tno+lUo+2ino+3) and 0(n o + 2,n o + 3) = 

I ln ln +iTno+2Tno+3)- Such eigenfunctions and their 
correspondent eigenvalues are obtained from the diago- 
nalization of the following tridiagonal sub matrix 

£(0) + r + s r 
r +2r r 

r + r + 

Above r = B/(4A), s = B 2 /[A{BA+g)\ and^W = 2£ x + 
BA + 2g. The difference in energy between the diagonal 
element in the middle of the matrix and the diagonal 
elements at the edges, B/(4A)-B 2 /{4A(BA + g)}, exists 
because state </>(no, no + 1) and 0(no + 2, no + 3) make 
virtual transitions to states that are out of the 'sub' chain 
created by the defects. 

The eigenvalues are therefore 



Ex = E<® + 



E, = E«» 



E, = E^ 



AB 2 A + 3Bg - Bu 
8A{BA + g) 

B(2BA + g) 

AA(BA + g) 

4B 2 A + 3Bg + Bu 
8A{BA + g) 



where u = y / 8B 2 A 2 + \%BAg + 9g 2 . 

To generate a W state we have to prepare an initial 
state with excitations on no + 1 and no + 2, since this 
state is the only one of the three which is connected to 
the two others in second order. The probability in time 
to obtain the initial state </>(n + l,n + 2) is 



P. 



0(ra o + l,n o +2) 



(*) 



1 + cos [(Es-E^t] 



(10) 



The probability to find state 4>{ n 0i tiq + 1) or state 0(n o - 
2, no + 3) is the same and is given by 



'n ,n + l 



)(*) = 



1 - cos[(£ 3 - E x )t] 



ill) 



The W state appears at the following instants of time 
(see Fig.l) 



(-l) fc arccos(-l/3) + 2ir[k - int(fc/2)] 
E z -E 1 • 



(12) 
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FIG. 1: We considered a chain with 12 qubits, B = 1 and 
A = 10. The defects are located on sites no, no + 1, no + 2 
and no + 3, where no = 3 and g = BA. The solid line gives the 
probability in time to find the initial state 0(no + 1, no + 2). 
Both probabilities to find state 0(no, no + 1) and state <j>(no + 
2, no + 3) coincide and they are given by the dot-dashed line. 
The vertical dashed lines correspond to the instants of time, 
Tk = Btk, where we have a W state. 

where k is an integer and int(fc/2) correponds to the in- 
teger part of the ratio k/2. 

As pointed out in the previous section, for the maxi- 
mally entangled state to be kept this way, at the moment 
where it is generated, the level spacings of the qubits 
involved in the process should become different. This 
detuning should be larger than the interaction strength 
among them. 

The W state with the bound pairs <^>(no, n + l), (f>(riQ + 



I, no + 2) and (j>(no + 2,n + 3) can also be created with 
only two defects located on sites n — 1 and n + 4, but 
here the situation is more delicate. As before, g has to 
be larger than the bandwidth of the bound pair band, 
but it cannot be close to BA, because this would create 
resonances with states that have one excitation on the 
defect [13. 



III. CONCLUSION 

We have investigated how the defects of a spin chain 
with strongly anisotropic coupling can be used to ob- 
tain EPR and W states. These are the states used in the 
study of bipartite and threepartite entanglement, respec- 
tively. We considered the XXZ model with defects, for 
this is the model used to describe some quantum com- 
puters, such as the one based on electrons on helium 
It was shown that even though the interaction among 
qubits is on all the time, by controlling the level spacings 
we can determine among which qubits the interaction is 
actually effective. This allows entangling only certain 
chosen sites. 

Over the years, magneto-chemists have refined the 
art of designing and growing crystals of quasi-one- 
dimensional magnetic materials. These systems should 
therefore be useful in the study of entanglement. 
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